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The  free  electron  laser  instability  is  investigated  for  a 

*  i 

relativistic  annular  electron  beam  propagating  through  a  helical  wiggler 
magnetic  field.  It  is  assumed  that  the  beam  is  thin,  with  radial 

n 

thickness  (2a)  much  smaller  than  the  beam  radius  (R  ) ,  and  that 

l  0 

v/y.  «  1,  where  v  is  Budker’s  parameter.  The  stability  analysis  is 
carried  out  within  the  framework  of  the  linearized  Vlasov-flaxwell 
equations  for  perturbations  with  general  azimuthal  harmonic  number  H. 
and  radial  mode  number  s,  including  the  important  influence  of  (a) 
finite  beam  geometry  in  the  radial  direction,  (b)  positioning  of  the 

r  r 

beam  radius  relative  to  the  outer  conducting  wall  (Rq/Rc),  and  (c)  finite 
wiggler  araplitude^C^B) .  All  of  these  effects  are  shown  to  have  an 
important  influence  on  stability  behavior.  Moreover,  the  maximum 
coupling  between  the  transverse  and  longitudinal  modes  increases 
substantially  with  increasing  radial  mode  number  s.  It  is  also 
found  that  the  transverse  magnetic  (TM)  mode  has  slightly  larger 
growth  rate  than  the  transverse  electric  (TE)  mode. 


■4  7 

■A  -  / 


■r—T 

V* 


3 


unclassified _ 

SECURITY  CLASMEICATION  OF  Twit  PACK  Data  SllffO 

REPORT  DOCUMENTATION  PACE  sKroeVcoMPLETmcVoRH 

T."  i^oUT  mumHH  '  [a.  qovt  accession  wo.  *.  recipient's  catalog  number 

IH  SAI-254-80-626-LJ  f  ^ _ 

■FST  ESanEggS  X  m~  tvre  ftp  «g»QWT  *  wtmao  eovCRgo 


BPOWT  MUMfl 


*  SAI- 2  54-80^-6  26-LJpl 


>/  JREE  g^ECTRON  LASER  ^STABILITY  FO! 

i  Relativistic  ^nnular ^electron  beam 

A'  HELICAL  WIGGLER  FIELD  ,  v . . 


FOR  A 
AM  IN 


Han  S . / Uhm 
Ronald  C.l Davidson 


I  Technical  Hep^t .  [ 

X-  I  - L  - 

RCRFORMINO  OSS.  RENORT  NUMBER 

PRI-17  ^ 

a.  CON  A  ACT  OR  CRAnV  NUMptNC*; 

>^ONR  ^^14-79^-^555/^ 


T  piiroSSo  oroanixation  namk  ano  aoorcss  "  TA$* 

ARIA  A  WORK  UNIT  NUMBERS 

Science  Applications ,  Inc. 

Plasma  Research  Institute 

934  Pearl  Street.  Boulder.  Colorado  80302 -  —  .-j - - - 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  r*-  REWIRT  DATE  ■■  J - - 

Office  of  Naval  Research  t  I  J  / 

Physics  Program  Office  r *Bt*  ’ 

14  MONlTO  ^IhVaGCNCY  ^N  AMS  •^AOORCIVU  3J5H57  SSS  CmSSSm  OHIcb)  IS.  SECURITY  CLASS.  (•!  thlm  r BBBli) 


unclassified 


Sa.  OECL  assi  ri  C  ATION/ OONNCN  AOINO 
SCHEDULE 


I  IS.  OISTRiaUTION  STATEMENT  f of  Nila  «lfM| 


Approved  for  public  release:  distribution  unlimited 


1 17.  DISTRIBUTION  STATEMENT  (•!  *•  B*BlrBtt  mnlBrB*  In  Black  M.  It  mttrn rail  Hmm  Rm^) 


I  is.  supplementary  notes 


IS.  KEY  BOROS  (CBmttm m  «  nmn  If  mmmnf  «4  M MrtJV  *r  M««*  nmN*I 

relativistic  electron  beam  for  free  electron  laser 


M.  ANSI  RACT  (CmmthmB  m  rBBBTBB  i 

The  free  electron  lase 
electron  beam  propaqat 
assumed  that  the  beam 
beam  radius  (Ro),  and 
stability  analysis  is 
Maxwell  equations  for 
and  radial  mode  number 


r  instability  is  investigated  for  a  relativistic  annular 
inq  through  a  helical  wigqler  magnetic  field.  It  is 
is  thin, with  radial  thickness  (2a)  much  smaller  than  the 
that  v/Yb  1,  where  v  is  Budker's  parameter.  The 
carried  out  within  the  framework  of  the  linearized  Vlasov 
perturbations  with  general  azimuthal  harmonic  number  l 
s,  including  the  important  influence  of  (a)  finite  beam 


COITION  OR  I  NOV  SS  IS  I 

S/M  0ISMJ-0I4-  SSOI 


unclassified 

SECURITY  CLASSIFICATION  OR  THIS  RR  (I 


SECURITY  CLASSIFICATION  OF  TNIS  FAST  (Dim  Data  BaMr*< 


20.  Abstract,  continued. 


geometry  in  the  radial  direction,  (b)  positioning  of  the  beam  radius  relative 
to  the  outer  conducting  wall  (R  /R  )  and  (c)  finite  wiggler  amplitude  (SB) . 
^11  of  these  effects  are  shown  So  Save  an  important  influence  on  stability 
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I.  INTRODUCTION 

There  is  a  growing  literature  on  the  free  electron  laser*  ^ 

which  generates  coherent  electromagnetic  radiation  using  a  relativistic 

electron  beam.  For  the  most  part,  previous  theoretical  analyses  of  this 

instability  have  been  carried  out  for  an  electron  beam  with  uniform 

3“6 

density,  neglecting  the  influence  of  finite  radial  geometry. 

Strictly  speaking,  a  more  accurate  theoretical  model  of  radiation 
generation  by  the  free  electron  laser  instability,  including  a  determina¬ 
tion  of  the  optimum  value  of  the  beam  radius  Rq,  requires  a  linear 
stability  analysis  for  perturbations  about  an  annular  electron  beam 
propagating  in  a  helical  wiggler  magnetic  field.  This  paper  develops 
a  self-consistent  theory  of  the  free  electron  laser  instability 
for  a  relativistic  annular  electron  beam  propagating  in  a  helical 
wiggler  field,  allowing  for  perturbations  with  general  azimuthal 
harmonic  number  i.  The  present  work  extends  the  previous  analysis^ 
by  the  authors,  carried  out  for  an  undulator  (multiple  mirror)  wiggler 
field  and  l  =  0  perturbations,  to  the  case  of  a  helical  wiggler  field 
and  arbitrary  azimuthal  harmonic  number  g. 

The  analysis  is  carried  out  within  the  framework  of  the  Vlasov- 
Maxwell  equations  for  a  relativistic  annular  electron  beam  propagating 
in  the  combined  transverse  wiggler  and  uniform  axial  guide  fields 
described  by 

3  B.e  -  6Bcosknz6  -  6Bsink«ze  , 

%  (Kz  0  %x  -  0  'vy 

where  Bq  and  6B  are  constants,  Aq  •  2n/kQ  *s  ax*8^  wavelen8th 

of  the  helical  wiggler  field,  and  6,6,  and  e  are  unit  vectors 

'vx  -vy’  -\,z 

in  the  x-,  y-,  and  z-d i rect i ons ,  respectively.  It  is  assumed  that 
the  beam  thickness  (2a)  is  much  smaller  than  the  mean  beam  radius  Rq, 
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and  that  v/y^  «  1»  where  v  is  Budker's  parameter,  and  y^mc  is  the 
characteristic  energy  of  the  electron  beam.  In  Sec.  II,  equilibrium 
properties  are  calculated  for  the  choice  of  equilibrium  distribution 
function  [Eq.  (14)] 


fj|  -  K  x  6(Cb  -  Cq) 6 (Cx  -  y^2Vq)G(Cz)  , 

Q 

where  C,  ,  Cx,  and  C  are  the  helical,  transverse,  and  axial  invariants 
n  2 

defined  in  Eqs.  (12),  (11),  and  (13),  respectively,  K  is  a  normalization 

constant,  and  Vq  and  Cq  are  constants  related  to  the  radial  thickness 

and  mean  radius,  respectively,  of  the  beam. 

The  formal  stability  analysis  for  the  free  electron  laser  instability 

is  carried  out  in  Sec.  Ill,  and  a  complete  dispersion  relation  [Eq.  (68)] 

is  obtained,  assuming  that  In)  -  (k  +  nk„)g,  cl  <<  k,.B,  c,  w  ,  where  u  and 

L;  D  u  d  c 

k  +  nkQ  are  the  complex  oscillation  frequency  and  the  wavenumber  of  the 

perturbation,  respectively,  =  cBq/y^c  is  the  electron  cyclotron 

2  1/2 

frequency,  c  is  the  speed  of  light  in  vacuo,  and  B^c  =  (yb  -  1)  c/y^ 
is  the  mean  axial  velocity  of  the  electron  beam.  The  resulting 
dispersion  relation  [Eq.  (68)]  constitutes  one  of  the  principal 
results  of  this  paper  and  can  be  used  to  investigate  stability  properties 
for  a  broad  range  of  system  parameters. 

In  Sec.  IV,  the  dispersion  relation  for  the  free  electron  laser 
instability  is  investigated  for  u  and  k  +  nkQ  near  the  simultaneous 
zeroes  of  the  transverse  (vacuum  waveguide)  dispersion  relation. 


*2  ‘  (k  +  nk0  "  V  jRc 

c  ; 


£«.. . 
Ki.s . 


TE  mode, 
TM  mode. 


and  the  longitudinal  dispersion  relation  r  +  (7r/2v,  )y  =  0, 

l ,n  ' b  *n,n 

where  R  is  the  radius  of  the  conducting  wall,  r  and  v  ,  are 
c  t  ,n  *n,n 

defined  in  Eqs.  (54)  and  (81)  and  B. ,  and  a  ,  are  the  s'th  roots 

x  , s  i  ,s 
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of  the  Bessel  function  J  , (6  ,  )  *  0  and  its  derivative  J',(o  ,  )  «  0, 

X.  1  ,  S 

respectively.  The  abbreviations  TE  and  TM  refer  to  transverse  electric 
and  transverse  magnetic  polarizations,  respectively.  It  is  shown  that 
the  maximum  coupling  between  the  longitudinal  and  transverse  modes 
occurs  at  a  value  of  Rq  satisfying 

Ro/R  /aM/a««.s  •  IE"°de> 

0  c  1  "°de  • 

Moreover,  the  coupling  coefficient  increases  considerably  with  increasing 
radial  mode  number  s.  It  is  also  found  that  the  TM  mode  is  slightly 
more  unstable  than  the  TE  mode. 
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II.  EQUILIBRIUM  CONFIGURATION  AND  ASSUliPTIONS 

The  equilibrium  configuration  consists  of  a  relativistic  annular 
electron  beam  propagating  in  the  combined  transverse  wiggler  and  uniform 
axial  guide  fields  described  by 


*  =  Bolz 


+  6B 


=  B.e  -  6Bcosk„ze  -  6Bsink„ze  ,  (1) 

(Kz  0  'voc  0 

where  B^  and  6B  are  constants,  \q  =  2tt /k^  =  const,  is  the  axial 
wavelength  of  the  helical  wiggler  field,  and  e^,  (^  ,  and  are 
unit  vectors  in  the  x-,  y-,  and  z-directions ,  respectively.  The 
electron  beam  has  a  characteristic  radial  thickness  2a  and  mean  radius 
Rq,  and  is  located  inside  a  grounded  cylindrical  conducting  wall 
with  radius  Rc>  We  introduce  cylindrical  polar  coordinates  (r,0,z) 
with  z-axis  along  the  axis  of  the  beam;  r  is  the  radial 
distance  from  the  z-axis,  and  6  is  the  polar  angle  in  a  plane 
perpendicular  to  the  z-axis.  In  cylindrical  coordinates,  the  trans¬ 
verse  wiggler  field  6B  in  Eq.  (1)  can  be  expressed  as 


6B 

% 


6B  e  + 
ro-r 


6B  e 
0^0 


=  -6Bcos(0  -  kQz)£r  +  SBsin (0  -  kQz)£6  , 


(2) 


where  ^  and  are  unit  vectors  in  the  r-  and  0-  directions,  respectively. 
To  make  the  analysis  tractable,  the  following  simplifying  assumptions 
are  made. 

(a)  The  thickness  of  the  annular  electron  beam  is  much  smaller 
than  its  mean  radius,  i.e., 


a/RQ  «  1. 


(3) 


6 


(b)  The  characteristic  transverse  momentum  of  a  beam  electron 
is  small  in  comparison  with  the  characteristic  directed  axial  momentum 
YfamVb,  i.e., 

k±\  «  Yb“Vb  £  pz  ,  (A) 

where  Vb  is  the  mean  axial  velocity  of  the  electron  beam. 

(c)  It  is  also  assumed  that 


v/Yfa  << 


(5) 


2  2  2 
where  Ybmc  is  the  characteristic  electron  energy,  and  v  =  Nbe  /me 

is  Budker's  parameter.  Here,  c  is  the  speed  of  light  in  vacuo,  -e  and 

m  are  the  electron  charge  and  rest  mass,  respectively,  and 

,2* 


N, 


de 


f  dr  r  n^  is  the  number  of  electrons  per  unit  axial  length. 

Jo  b 


Consistent  with  Eq.  (5),  the  equilibrium  self  fields  can  be  neglected 

in  comparison  with  the  applied  magnetic  field  in  Eq.  (1). 

(d)  It  is  assumed  that  the  axial  wavelength  \ ^  =  2ir/kQ  of  the  helical 

wiggler  field  is  sufficiently  short  that 

2 


l“0  '  “c 

where  and  u>  are  defined  by 
0  c 


61T  ,2  2 

^X'‘V 


(6) 


eBr 


k_.V,  and  u  =  - 

Ob  c  YKmc 


(7) 


Consistent  with  Eq.  (6),  it  is  also  assumed  that 

* 


2u  oj. 

c  0 


^bmVb 


6B 

n, 

bT 


«  1  , 


(8) 


i  ”0  c 1 

which  is  easy  to  satisfy  for  sufficiently  short  wiggler  wavelengths. 


(e)  In  the  stability  analysis,  it  is  assumed  that  the  wave 
perturbations  are  far  removed  from  resonance  with  the  transverse 
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cyclotron  motion  and  the  axial  motion  of  the  beam,  i.e., 


1“  '  <k  +  nk0)VJ  <<  “c,k0  \  * 


(9) 


where  w  is  the  complex  oscillation  frequency,  and  k  +  nkg  is  the 
wavenumber  of  the  perturbation.  It  is  also  assumed  that  the  cyclotron 
motion  and  the  axial  motion  is  nonresonant  with 


,2  2  ,  2 

k0vz  4  “c 


(10) 


where  =  pz/ym  is  the  axial  velocity  of  a  typical  beam  electron. 

Within  the  context  of  Assumption  (c),  there  are  three  exact 

g 

invariants  associated  with  the  single-particle  motion  in  the  equilibrium 
field  JgG  [Eq.  (1)].  These  are  the  transverse  invariant  C^, 


2  .  2  ,  2eB0  ,  „  ,  .  2e 

6  +  (pz  "  V"V  +  d^'6S’ 


zi  =  Pr  +  Po  + 


(ID 


the  helical  invariant  C^, 


Ch  =  P0  +  ^  (PZ  ~  VV  +  ^  rsin(e  "  koz)  * 


(12) 


and  the  axial  invariant  C  defined  by 

z  * 


(13) 


where  ^  =  (p^,  p0,  pg)  =  ymv  is  the  mechanical  momentum,  P0  =  r(pfl  -  eBQr/2c) 

is  the  canonical  angular  momentum  associated  with  the  axial  field  Bq, 

2  24  22  1/2 

ymc  =  (m  c  +  c  £  )  is  the  relativistic  electron  energy,  and  y^  = 

(1  -  V2/c2)"1/2  =  const. 

D 

In  the  present  analysis,  we  consider  the  general  class  of  annular 

g 

electron  beam  equilibria  described  by 

fb  =  K  6<Cx  -  ^bm2vO)6<Ch  '  C0)G(Cz)  * 


where  v^  is  a  constant  related  to  the  radial  thickness  of  the  beam, 


(14) 
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defined  by 


9 


In  the  lowest  order  calculation  presented  here,  the  electron  density 


profile  in  Eq.  (19)  is  independent  of  6  and  z  for  system  parameters 
satisfying  Eqs.  (3)  and  (6).  Additional  equilibrium  properties  assoc¬ 
iated  with  the  distribution  function  in  Eq.  (14)  are  discussed  in  Ref.  8. 
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III.  LINEARIZED  VLASOV -MAXWELL  EQUATIONS 

In  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwell 
equations  to  obtain  an  eigenvalue  equation  describing  the  free  electron 
laser  instability  in  an  annular  electron  beam.  In  the  subsequent 
analysis,  it  is  assumed  that  all  perturbed  quantities  have  temporal 
and  spatial  variations  of  the  form 


6i|)(x,t)  =  6tKx)exp(-iut)  , 


where  Imui  >  0,  and  ,p(x)  is  the  amplitude  of  the  perturbation.  Using 
the  method  of  characteristics,  and  neglecting  initial  perturbations, 
the  perturbed  electron  distribution  function  is  given  by 

ft 


ifb(x^,t)  =  e 


dt’ 


1  3 


-v 1 6<j>(x' , t * )  -trr  W’f) 


v'  X  V'  X  6A(x' , t' ) 


(21) 


•^bV’^ 


where  dx'/dt'  =  v'  and  d^'/dt'  =  -ev'  x  B  (x')/c,  and  the  particle 
trajectories  (x' ,  ^')  in  the  equilibrium  field  configuration  satisfy 
the  "initial"  condition  x'(t'  =  t)  =  and  ^’(t'  =  t)  =  £.  In  Eq.  (21), 
the  perturbed  electric  and  magnetic  potentials,  64  and  6A,  are  determined 
self-consistently  from  the  linearized  Maxwell  equations 


(•'■?&) 


641  =  Ane  I  d  p  6f^  , 


(22) 


and 


(■■ -?&)**•  “I 


d  p  v  6f  , 

^  D 


(23) 


respectively,  where  use  has  been  made  of  the  Lorentz  gauge  condition 


1  3 


V  •  M(x,t)  +  -  --  H(x,t)  =  0  . 


V* 


1 


V 


(24) 
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-  [  (<Wk0)sin(e  -  k0z)rosu;cT  -  vzcos(e  -  k0z)sinW(;T]}  ,(28) 

where  use  has  been  made  of  Eqs.  (12)  and  (26),  and  t  =  t'  -  t. 

As  evident  from  Eqs.  (27)  and  (28),  the  radial  and  azimuthal 
orbits  contain  oscillatory  contributions  proportional  to  cosh^t, 
simo^t,  sin(e  -  k^z  -  k^v^)  and  cos(6  -  k^z  -  ^qv2t)«  For  present 


12 


purposes,  in  the  t 1 -integration  on  the  right-hand  side  of  Eq.  (21), 
we  retain  term?  proportional  to  [to  -  (k  +  nkQ)pz/Ym]  \  and  assume  that 
the  value  of  [w  -  (k  +  nkQ)pz/Ym]  is  well  removed  from  resonance 
with  the  cyclotron  motion  and  the  axial  motion  of  the  beam  [Eq.  (9)], 


|w  -  (k  +  nk0)pz/Ym|  «  wc  ,  kQVb  . 

Within  the  context  of  Eqs.  (1)  and  (9),  the  perturbed  distribution  function 
in  Eq.  (21)  can  be  approximated  by^ 


fb(x,£)  =  -e  dTexp(-iu)i)  \2  I  ymitu ( 6 4>  -  —  •  6£) 


•  k  4)M ^  -  h'  •n'tj-ir 

1  d  pj_  V7. 

2  222  1/2 

where  y  =  H/mc  =(l+^/mc)  .  Moreover,  making  use  of  Eq .  (9), 

only  those  contributions  to  v^  and  v^  proportional  to  sin( 0  -  k^z  -  ^qvzt) 

and  cos (6  -  k  z  -  kQv  r)  are  retained.  That  is,  on  the  right-hand 
U  z 

side  of  Eq.  (29),  we  retain  contributions  to  v^  and  v^  of  the  form 

v<  =  v  - S. -  i!  cos(6  -  k  z  -  k.v  t )  ,  (30) 

r  z  —  u)  B_  0  U  z 

U  c  0 

and 

v>  =  _v - 9 - sin(e  -  knz  -  k.v  t)  .  (31) 

0  z  Wq  -  o)c  BQ  0  0  z 

Finally,  since  the  oscillatory  modulation  of  the  radial  and  azimuthal 
orbits  is  small-amplitude  [Eqs.  (3),  (6),  and  (9)),  we  approximate 


r’  =  r  ,  e’  “  6  , 


in  the  arguments  of  the  perturbation  amplitudes  on  the  right-hand  side 


of  Eq.  (29),  i.e., 


i(^n^  (r '  )oxp(Ue'  )  ~  (r)exp(ite) 
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Substituting  Eqs.  (25)  and  (30)  -  (32)  into  Eq.  (29),  and  making 
use  of  Eq.  (9),  we  find  after  some  straightforward  algebra  that  the 
perturbed  distribution  function  6f^(^t,^)  can  be  approximated  by 


6fb(^*^  =  e  l 


exp {i f  1.6  +  (k  +  nkQ) z] }  j  /  (r)\ 

,n  oi  -  (k  +  nkQ)vz  |n^l  c  z,l) 


+  1  1*  -_-c- _ 11  [\  A(n+1)  -  X  A(n-1)  1 

2  c  oj„  -  Bq  ln+lA6,i-l  n-1  6,1+1 J 


+  \  —  — ^ - -  I5-  [X  A 

2  c  a)  -  B-.  n+1  ' 

c  0  0 


+  x  lA<n;i>]V  • 

n,£-l  n-1  n,£+l  / 


where  the  function  X  ,  (]j,a),k)  is  defined  by 

8fb  sfb 

Xn,(^,u,k)  =  2trmw  -  (k  +  n’k0)Pz]  — j  +  (k  +  nkQ)  — -  . 

It  is  convenient  to  introduce  the  dimensionless  potential 
amplitudes  defined  by 


(33) 


(34) 


Vr>  -^2  ^n)w  -r  Az°lM!  ■ 

me 


(35) 


;  /_\  e  (n+1) ,  , 

ve,i±i(  ^  2  Ae,i±i(  ^  ’ 

me 


and 


—  A(n+1)(r) 

2  r , !±1 'r'  ‘ 
me 


Ar,l±l(r)  = 

Moreover,  we  introduce  the  dimensionless  parameter 


(36) 


(37) 


A  = 


e6B 


2^bmC  k0 

and  define  the  effective  susceptibility  xn  n<  (u>,k)  by 
,  ,  ,  ,  2fRc  .  f  ,3  Xn,(S’u,k) 

Vn’(u’k)  =  4ne  j0  rdT  p  • 


(38) 


(39) 
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To  simplify  the  present  analysis,  it  is  also  assumed  that 


w0  “c  * 


(40) 


where  ^  kgV^  and  wc  ”  e®0^Ybmc’  EQuatl°n  (40)  is  easily  satisfied 
in  parameter  regimes  of  practical  Interest  for  free  electron  laser 
applications. 

f  3 

The  integral  jd  pv6f^  in  Eq.  (23)  is  required  to  evaluate 
the  r-,  0-,  and  z-components  of  the  perturbed  current  density. 

Within  the  context  of  Eqs.  (3),  (6),  and  (9),  It  is  important  to 
note  from  Eqs.  (11)  and  (14)  that  the  quantity  A  defined 

g 

in  Eq.  (34)  is  an  even  function  of  p^  -  PfQ.  where 

U) 

p  n  »  2>  me  - - -  Acos  (0  -  k  z)  .  (41) 

ru  b  -  a)  0 

0  c 

Therefore,  it  follows  from  Eq .  (33)  that 


.3 

d  p  v 


d  P(Pr0/>bm)Afb 


where  we  have  approximated  >  -  which  is  consistent  with  Eq.  (9). 

That  is,  the  perturbed  radial  current  is  given  approximately  by 


f  .3 

-e  d  p  v 


-2ecAcos(t’  -  k,.zKn. 

l)  b 


(42) 


where  use  has  been  made  of  Eq.  (4ni.  It  can  also  be  shown  that  the 
perturbed  azimuthal  current  is  given  approximately  by 


-e 


2eo\sin(e  -  kgzJfnj,  . 


(43) 


Making  use  of  Eqs.  (33)  -  (43),  the  Maxwell  equations  (22)  and 
(23)  can  be  expressed  in  the  approximate  form 


v 


« 


v 
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/,  ,  3  ,2  2  2\  .  «(r  -  RQ) 

Srbrf7'A2+^  -  <k  *  nk0)  /♦t(r)  ’  - 2 -  °(“*k)  *  (44) 

l  r  C  /  YbR0 

{1  3  3  U  +  l)2  +  1  .  w2  .  .  ...  ,2\  t 

r  3  r  r  3 r  -  [k  +  (n  -  l)kQ]  /Ar> 


_  2H}±11  a  (r)  .  1( ,r "  r.q) 

2  e,£+iv  ;  r 


A  o  (uj,k) 


3  a  («+ir  +  i  ^  ,,  ,  ,  , ,,  ,2V  ;  , 

3 r  r  3 r  r 2  £2  "  *k  +  (n  1)k0^  /Ae,t+l(r) 

*  -  2  ‘r.-lW  ‘  -L-(I‘I'WI  ’ 

r  0 

i7  r  r,7  ^2  +  “2  -  Ik  +  (n  +  l)k0]  \Ar>  c_i(r) 


.  2J_Ci!U  A  {r) , 6 (r v, 

2  A0,t-lX '  R. 


Ao(ui.k) 


y  , r  r  .•  7  *  ~  S  +  c2  lk  +  (n+1)k0]  J  Ae,e-i(r) 

"Mfi-l)*  4<r  -  ty 

+  T  -  \,t-l(r)  ■  -  Rq—  ('iA)c’(“-k)  •  (4t 

for  a  thin  annular  beam  satisfying  Eq .  (3).  In  Eqs.  (44)  -  (48),  the 
sourre  function  i'(u,k)  is  defined  by 


(“,k)  =  xn,nW  L^txn,n+lAc,i-l(R0)  '  xn  ,n- 1%  ,  t+1  (R0)  1 


"  A  lxn.n+lAr,C-l(R0)  +  xn  ,n- lAr ,  t  +  1 ( R0)  1  ' 

Since  the  right-hand  side  of  Eq .  (44)  vanishes  except  at  r 
Eq.  (44)  can  he  expressed  as 


(\  3  3  82  2\ 

\r  3r  r  r  ~  2  pn  J 


4t(r)  3  0  » 
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for  r  i  Rg.  In  Eq.  (57),  the  parameter  pn  is  defined  by 

2  2,2  „  ^  .  .2 

Pn  =  w  /c  -  (k  +  nk0)  . 


The  solution  to  Eq.  (50)  is  given  by 


*„<r) 


J  AVpi>t)  +  “V* nr)  •  R0<t-<Rc' 

l  CJi(pnr)  •  0  r  <  R0  * 


where  J^(x)  and  N  (x)  are  Bessel  functions  of  the  first  and  second 
kind,  respectively,  of  order  t,  R^  is  the  radius  of  the  conducting  wall, 
and  the  constants  A,  B,  and  C  are  determined  from  the  appropriate 
boundary  conditions.  Multiplying  Eq .  (44)  by  r  and  integrating 
from  Rq(1-c  )  to  Rq(Ht),  with  e  ■*  0+,  we  obtain^ 


Ft  ,n(u,k)4,«.^R0)  “  2  a  (a)  ,k)  , 


where  the  longitudinal  dielectric  function  ^(iD.k)  is  defined  by 

Jp  <C  )/J.  (£„) 

r£*"(W,k)  =  WW  ~WW  ’  ^ 

with  f,  =  p  R„  and  t,  =  R  p  . 
n  *n  0  n  crn 

The  transverse  dielectric  functions  associated  with  Eqs.  (45)  -  (48) 
are  calculated  in  detail  in  Appendix  A.  For  example,  for  0  <  r  <  Rq, 
the  solutions  to  Eqs.  (45)  and  (46)  are  given  by 

\,m<r)  •  c"+1)jt(P„-ir>  +  c!‘+UWVlr>  • 


v«+i<r>  ■  -  icri>jl+2<p„-ir)  • 

2  2  2  2  (t+1) 

where  p  ,  =  u)  /c  -  [k  +  (n-l)/k-J  ,  and  the  constants  C'  are 

n- 1  u  ± 

determined  from  the  appropriate  boundary  conditions.  After  some 
st raight forward  algebraic  manipulation,  which  is  summarized  in 


v 
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Appendix  A,  it  can  be  shown  from  Eqs.  (45)  and  (46)  that 


Dt+l(“,k)C-1+1>Jl(Cn-l)  “  IAo<“,k)  ’  (5' 

where  the  transverse  dielectric  functions  D++^(w,k)  and  D^+^(w,k)  are 
defined  by 


VlJI+l(Cn-l)Jhl(tn-l) 

VWWW  ’ 


8E+l(;n-l)/l,J£^n-P 


Here,  the  prime  (')  denotes  (d/dx) (x) .  In  Eqs.  (58)  and  (59), 
42,  =  [u,2/c2  -  (k  +  n'k0)2]R2,  £2,  =  [w2/c2  -  (k  +  n’k0)2]R2,  and 
the  function  , ( c)  is  defined  by 

2J£, (C)J’ , (C) 

8S-'(°  =  Jat(C)N^,(C)  +  ‘Ti'(t)Nt,(0  ’  1 

Similarly,  the  solutions  to  the  coupled  differential  equations 
(47)  and  (48)  are  given  by 


Vi-i(r>  ■  ““'“-Ww*  -  • 


for  0  <  r  <  Rq.  From  Eqs.  (47)  and  (48),  we  obtain 


D^1(w,k)C^£"1)Ja(Cn+1)  -jAoU.k)  , 


D;_1(„,k)c<‘-1)Jt_2(tntl) 


Ao(uivk)  , 


V* 
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where  the  transverse  dielectric  functions  D^^w.k)  and  D~  ^w.k)  are 
defined  by 


Dj^Cw.k) 


8l-l(Cn+l)/,t‘V^n+l* 

8£-l(i’n-»-l)Nl(('n+l) 


VW 


(64) 


and 


D^_l  (>J»k)  = 


;n-HJ£-l  ^n-H^  Jl-1  ^n+1^ 
Jl(f,n+l^JJt-2(^n+l^ 


Substituting  Eqs.  (55)  and  (61)  into  Eq .  (49),  the  source 
function  a(oj,k)  can  be  expressed  as 


(65) 


o<u,k)  -  V„VV  - 


*  ‘n,n+l'-;i'1>J«<Cn+l)1  • 


(66) 


Equations  (53),  (56),  (571,  (a*.),  and  (63)  constitute  a  set  of  linear 

*  (i-1) 

algebraic  equations  relating  ^^(Rq)»  C*  and  C+  .  However, 

we  note  from  Eq.  (66)  that  the  common  source  term  o(u),k)  is  related  only 
to  the  constants  ^^(rq)>  ,  and  In  this  regard,  the 

equations  containing  and  ^  are  completely  decoupled  from 

tlie  remaining  equations.  Without  loss  of  generality,  Eqs.  (56) 
and  (63)  are  therefore  omitted  from  the  subsequent  stability 
analysis.  Making  use  of  Eqs.  (53),  (57),  and  (62),  we  obtain  the  matrix 
equation  relating  ^(Rq),  1^Jn^n+l^’  and  C- £+^  J£  ^n-l^  ’ 


f 


\ 
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which  is  similar  in  general  form  to  the  matrix  equation  obtained  by  Uhm 
and  Davidson^  for  the  free  electron  laser  instability  associated  with 
an  annular  electron  beam  propagating  through  an  undulator  multiple 
mirror  magnetic  field.  Of  course,  the  transverse  dielectric 
functions  , ,  and  D+  ,  for  a  helical  wiggler  field  are  different  from 
those  for  an  undulator  field. 

The  condition  for  a  nontrivial  solution  to  Eq.  (67)  is  that 
the  determinant  of  the  matrix  vanish.  This  gives  the  general  dispersion 
relation 


D2-lD2+l^r2,n  +  z_2  Xn,nj 


(68) 


"  -*2(D2-l*n,n-l  +  D2+lXn,n+l)r2,n  =  °  ’ 


which  can  be  used  to  determine  the  complex  eigenf requency  u  in  terms  of 
k  +  nkg,  k^,  v.  A,  and  A . 
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IV.  FREE  ELECTRON  LASER  STABILITY  PROPERTIES 


We  now  investigate  the  free  electron  laser  stability  properties 

predicted  by  Eq.  (68)  for  a  relativistic  annular  electron  beam 

propagating  in  a  helical  wiggler  field.  Consistent  with  Eqs.  (5) 

and  (9),  the  eigenf requency  tu  can  be  approximated  by  w  -  (k+nkg)^. 

22  2  1/2 

therefore  approximate  G  ,,  =  {<*»  /c  -  [k+(n+l)krt]  }  R  by 

n+i  u  c 


We 


^n+1  1 


t  (k+nk  )Z 

- 2 -  +  2k0(k+nkQ)  +  kQ 


1/2 


(69) 


After  careful  examination  of  Eq .  (64)  together  with  Eqs.  (60)  and  (69), 

it  is  evident  that  the  transverse  dielectric  function  D*_^(u,k)  is 

never  equal  to  zero  for  k  +  nk^  >  0.  In  this  regard,  for  small  wiggler 

amplitude  (A  «  1) ,  we  investigate  free  electron  laser  stability 

properties  for  w  and  (k  +  nk^)  near  the  simultaneous  zeros  of  the  (A  •+  0) 

transverse  dispersion  relation  D^+^(u,k)  =  0,  and  the  (A  -*■  0)  longitudinal 

2 

dispersion  relation  r  +  (ir/2v.  )x  =  0.  For  small  but  finite  A  , 
y  f,n  Tb,An,n 

the  general  dispersion  relation  in  Eq.  (68)  can  then  be  approximated 
by  the  simplified  form 


Dm(ri,n  +  *n,n^ 


=  -A 


*n,n-l^£,n 


(70) 


Equation  (70)  is  one  of  the  principal  results  of  this  paper  and  can  be 
used  to  investigate  stability  properties  for  a  broad  range  of  system 
parameters . 

Making  use  of 


pn  »  iqn  ;  i(k  +  nk0)/yb  , 


(71) 


the  longitudinal  dielectric  function  ^  can  be  expressed  as 

r,,„  ■  2  h<<*„V  2  i‘tr<.nVK(T,;-R0)  -  l>nVKARc>  ' 


(72) 


V‘ 
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for  k  +  nkg  >  0.  In  Eq.  (72),  I^(x)  and  K^(x)  are  modified  Bessel 
functions  of  the  first  and  second  kind,  respectively,  of  order  l,  and 
h(qnRo)  is  the  effective  longitudinal  wave  admittance.  For  short 
axial  wavelengths  with  (R^  -  Rq)  (k  +  nk^)  >>  Yb»  it;  is  readily  shown 
from  Eq.  (72)  that  the  wave  admittance  h  can  be  approximated  by 

k  +  nk~ 


h(«nV  =  2  -V 


R0  ' 


(73) 


For  free  electron  laser  applications  characterized  by  axial  wavenumber 
2 

k  +  nkg  =  (1  +  Eq.  (73)  constitutes  an  excellent  approximation. 

For  very  small  wiggler  amplitude  (A  -*■  0)  ,  or  for  very  low  beam 
density  (v  -*-0),  it  is  evident  from  Eq.  (70)  that  the  linear  dispersion 
relation  for  transverse  perturbations  is  given  by 


Vl(a)’k)  =  °’ 


(74) 


and  the  linear  dispersion  relation  for  longitudinal  perturbations  is 
given  by 


r£,n(w’k)  +r2  *n,n  =  ° 


(75) 


2Y 


where  k  +  nk^  >  0  has  been  assumed.  Making  use  of  Eqs.  (59)  and  (60), 
the  dispersion  relation  in  Eq.  (74)  can  be  expressed  in  the  equivalent 
form 


S  ~  (k  +  nk0  “  V' 

c 


Vl,s 


(76) 


for  the  transverse  electric  (TE)  polarization,  and 

2  62 
^2  ~  (k  +  nk0  "  ko^  =  ‘  » 

c  R 

c 


(77) 


for  the  transverse  magnetic  (TM)  polarization.  In  Eqs.  (76)  and  (77), 

a  ,  is  the  s'th  root  of  J',(a  ,  )  =  0,  and  6  is  the  s'th  root 

i  ,s  l  l  , s  t',s 
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of  J  , (8  ,  )  =  0.  Equations  (76)  and  (77)  are  the  familiar  TE  and  TM 

A#  X.  }  S 

mode  dispersion  relations  for  the  vacuum  waveguide  modes. 

For  future  reference,  Taylor  expanding  the  transverse  dielectric 
function  Di+^(w,k)  in  Eq.  (59)  about  the  vacuum  waveguide  results,  it  is 
straightforward  to  show  that  D^+^(m,k)  can  be  approximated  by 

-  ,  ,  v  RC  f  2  ,..,'2l  JJH-l(ttft+l,s) 

ii  (t»>,k)  “  ^2+1  s  ( £+1)  t  /„  R  /R  ) 

Z+l  2  2  £+i»s  i  *V  A+l.sW 

1+1,3  (78) 

'  i  ' 

uk  ,  ,  \2  ail+l,s 

x  ~2  ~  (k  +  nk0  -  k0)  -  “1 -  * 

c  R 

v  c 

for  the  TE  modes,  and 

-  ,  <  fJhi(Vi.s>  f 

,+i  ’  ' 2  vvi.sW 


(k  +  nk  -  k  )2  -  "^1>S 
R 

c 


for  the  TM  modes. 

In  order  to  investigate  the  influence  of  axial  momentum  spread 
on  the  free  electron  laser  instability,  we  assume  an  axial  distribution 
function  of  the  form 


G(C  )  =  - 


(C,  -  V*V  + 

where  A  is  the  characteristic  spread  in  about  the  mean  value  y^mV^. 
We  further  assume  that  the  spread  A  is  small  in  comparison  with 
y,  mV.  ,  i.e.,  A  «  y,mV  .  Substituting  Eqs.  (14),  (34),  and  (80) 

D  D  D  D 

into  Eq.  (39),  we  then  obtain 

„  co2  -  (k+nk  )(k+n’kn)c2 

2v  0  U  /0, 

X  ,  =  ‘  o  2  9 

n,n  Yb  [u  -  (k+nkQ)Vb  +  i|k+nk0|A/rbm] 
within  the  context  of  Eqs.  (6)  and  (9). 


4 
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For  small  but  finite  wiggler  amplitude  and  low  beam  density, 
we  make  use  of  Eqs.  (72),  (78),  (79),  and  (81),  to  simplify  the 


full  dispersion  relation  in  Eq.  (70).  This  gives 

2 


to^  /i  I  .  ,  v2  a£+l ,  s 

2  ^k  +  nk0  ~  k0^  2 


R 


2vc 

A 


k+nk„  A 


to  -  (k-Hnko)V^  +  i  — — - — 
Ybm 


2  £+1  q 

2k0(k+nk0)  -  k*  - 


R 


} 


(82) 


4A 


2  vc 


Y,k 

b  c 


ko(k+nko_ko)  ■ 


&+l,s 


R 


«ts(VRc>  • 


for  the  TE  mode  polarization,  and 

2  8^  l 

to  n  ,  ,  .  .  2  ft+l,s 

2  ~  (k+nk0  ko^  2 


2vc 

* 


2kQ(k+nkQ)  -  kQ 


to  -  (k+nkA)V  +  i 
U  b 


I  k+nk0  j  A 


Ybm 


m,s 


(83) 


=  4\ 


2  vc 


Vc 


6 


ko(k+nVko)- 

R 

c 


OVRc> 


for  the  TH  mode  polarization.  In  Eqs.  (82)  and  (83),  the  TE  and  TM 

coupling  coefficients,  (R.-/R  )  and  (Rn/R  )»  are  defined  by 
r  °  £s  0  c  J?.s  U  c 


a„ 


<4(VV  = 


am,s  "  (J+1)2  i  J*+i(a*+i,s)  J 


(84) 


and 


<s<VRc> 


fVW^sVVl 


Jt+liec+l,s^ 


(85) 


respectively.  Equations  (82)  and  (83)  are  the  dispersion  relations 
used  in  the  remainder  of  this  section,  and  can  be  used  to  investigate 
stability  properties  for  a  broad  range  of  system  parameters  of  experimental 
interest . 


i 


A  careful  examination  of  Eqs.  (64)  and  (85)  shows  that  the 
coupling  between  the  transverse  and  longitudinal  inodes  is  maximum 
whenever  J'(a  n  R../R  )  =  0  for  the  TE  mode,  and  •11(6.,,  Rft/R  )  m  0 
for  the  TM  mode.  Here,  the  prime  (')  denotes  dJ^(x)/dx.  In  this  context, 
it  is  found  that  the  maximum  growth  rate  for  perturbations  with 
azimuthal  harmonic  number  2  occurs  for  a  value  of  Rq/R given  by 


R0/Rc 


/a 

/e 


2+1 ,  s 
2+1,  s 


TE  mode  , 
TM  mode  , 


(86) 


where  a  is  the  first  root  of  J'(a  .)  =  0.  Equation  (86) 

X.  m  JL  A_  X*  y  X 

is  valid  only  when  a  ,  <  a.,,  for  the  TE  mode,  and  a„  ,  <8.,, 

J  2,1  -  2+1, s  2,1  -  2+1, s 

for  the  TM  mode.  For  a  >  a  (TE) ,  or  a  >  f5  (TM)  , 

X  y  X  X  I*  1.  j  S  A  ^  J.  X  I  X  y  S 

the  maximum  growth  rate  occurs  for  R„/R  =  1. 

0  c 

Shown  in  Fig.  1  are  plots  of  (a)  the  ratio  R^/R^  that  satisfies 

E  M 

Eq.  (86),  and  (b)  the  corresponding  coupling  coefficients  and 

for  several  values  of  azimuthal  and  radial  mode  numbers,  2  and  s. 

Several  points  are  noteworthy  in  Fig.  1.  First,  the  maximum  couplirj 

occurs  at  R../R  =  1  for  the  (2,s)  =  (0,1)  mode,  where  a„  .  =  3.83, 

u  c 

a  ~  1.84  and  p,  .  ~  3.83.  Second,  except  for  azimuthally  symmetric 

I  »  1  l  y  l 

perturbations  (i  =  0),  the  value  of  Rq/rc  corresponding  to  maximum 

coupling  increases  with  increasing  mode  numbers  2  and  s.  Third, 

from  Fig.  1(b),  for  a  specified  value  of  azimuthal  harmonic 

number  2,  the  maximum  cojpling  coefficients  corresponding  to  the  values 

of  Kq/rc  in  Eq .  (86),  increase  rapidly  with  increasing  radial  mode 

number  s .  Finally,  for  specified  values  of  (2,s),  we  note  from 

M  E 

Fig.  1(b)  that  the  maximum  value  of  Q‘  g  is  slightly  larger  than  g. 

I 


I 
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Defining  the  normalized  Doppler-shifted  eigenf requency  by 

ft  -  tu  “  (k+nk0)Vb]/k0c  ,  (87) 

we  calculate  the  normalized  growth  rate  =  lmft  numerically  from  Eqs.  (82) 

and  (83)  for  a  broad  range  of  system  parameters  y,  ,  v/y,  ,  A,  A/y  me, 

D  D  D 

Rq/R^,  £.,  and  s.  Shown  in  Fig.  2  are  plots  of  the  normalized  growth 

rate  fh  versus  (k  +  nk^/kp  [Eqs.  (82)  and  (83)]  for  (£.,s)  =  (1,3), 

y,  =  10,  v/y,  =  0.02,  A^  =  0.01,  A/y,  me  =  0.002,  and  Rn/R  =  a .  ./a  0 

b  b  b  Uci,iz,3 

for  the  TE  mode,  and  R^/R^  =  aj  3  ^or  t*ie  ™  i^ode.  It  is  evident 

from  Fig.  2  that  the  TM  mode  is  slightly  more  unstable  than  the  TE 
mode,  which  is  consistent  with  Fig.  1(b). 

Shown  in  Fig.  3  are  plots  of  the  normalized  TE  mode  growth  rate 
versus  (k  +  nk^/kg  obtained  from  Eq .  (82)  for  rq/rc  =  l^a£+l  s 
[Eq.  (86)].  The  two  plots  correspond  to  (a)  Z  =  1  and  several  values 
of  s,  and  (b)  s  =  3  and  several  values  of  £,  and  parameters 
otherwise  identical  to  Fig.  2.  As  predicted  in  Fig.  1(b),  it  is 
evident  from  Fig.  3(a)  that  for  a  specified  value  of  azimuthal  mode 
number  l,  the  growth  rate  and  range  of  k-space  corresponding  to  instability 
increase  rapidly  with  increasing  radial  mode  number  s.  For  R^/R^ 
satisfying  Eq.  (86),  we  therefore  conclude  that  perturbations  with 
high  radial  mode  numbers  exhibit  stronger  instability  than  the  fundamental 
mode  (s=l) .  Moreover,  the  k-value  corresponding  to  maximum  growth  is 
somewhat  reduced  as  the  radial  mode  number  s  is  increased.  For 
example,  in  Fig.  3(a),  the  maximum  growth  rate  occurs  at  k  +  nkQ  =  195  k^ 
for  s  =  1,  and  at  k  +  nkg  =  187  k^  for  s  =  5.  The  maximum  coupling 
inefficient  in  Fig.  1(b)  decreases  slowly  with  increasing  azimuthal 
mode  number  5,  as  does  the  instability  growth  rate  [Fig.  3(b)]. 

The  instability  results  for  the  TM  mode  are  similar  to  those  for  the 
TE  mode.  However,  the  TM  mode  is  somewhat  more  unstable  than  the  TE  mode. 
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V.  CONCLUSIONS 


In  this  paper,  we  have  investigated  the  free  electron  laser 

stability  properties  for  a  relativistic  annular  electron  beam  propagating 

in  combined  transverse  wiggler  and  uniform  axial  guide  fields.  The 

stability  analysis  has  been  carried  out  within  the  framework  of  the 

linearized  Vlasov-Maxwell  equations.  The  equilibrium  properties 

and  basic  assumptions  were  summarized  in  Sec.  II.  The  formal 

stability  analysis  was  carried  out  in  Sec.  Ill  for  general 

azimuthal  harmonic  number  £,  and  a  complete  dispersion  relation 

[Eqs.  (68)  or  (70)]  for  the  free  electron  laser  instability  was 

obtained,  including  the  important  influence  of  finite  radial  geometry. 

In  Sec.  IV,  this  dispersion  relation  was  investigated  for  u  and  k  +  nk^ 

near  the  simultaneous  zeros  of  the  transverse  and  longitudinal  dielectric 

functions.  It  was  shown  that  the  maximum  coupling  between  the  longitudinal 

and  transverse  modes  occurs  for  values  of  R_/R  satisfying  R„/R  = 

0  c  U  c 

a  ./a  .  for  the  TE  mode  polarization,  and  Rn/R  =  a  ./B  . 
for  the  TM  mode  polarization.  Moreover,  the  strength  of  the  coupling 
increases  considerably  with  increasing  radial  mode  number  s.  It  was 
also  found  that  the  TM  mode  is  slightly  more  unstable  than  the  TE  mode. 
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APPENDIX  A 


ELECTROMAGNETIC  POTENTIALS  IN  A  CYLINDRICAL  WAVEGUIDE 


In  this  appendix,  we  investigate  properties  of  the  electromagnetic 

potentials  in  a  cylindrical  waveguide  with  radius  R  .  Within  the  context 

c 

of  the  Lorentz  gauge  condition. 


V 


+ 


1  3__ 
c  3  t 


<Kx,t) 


0  , 


(A.  1) 


the  electromagnetic  potentials,  $  and  A,  are  determined  from  the 
Maxwell  equations 


t)  =  -4np(x,t)  , 


and 


('■-7  &)“<■ 

(v2  ~  7  =  ’  ^  ’ 


(A. 2) 


where  p(x,t)  and  J(x,t)  are  the  charge  and  current  densities, 
respectively,  and  c  is  the  speed  of  light  in  vacuo.  In  the  present 
analysis,  we  assume  that  all  quantities  vary  with  space  and  time 
according  to 


i|/(x,t)  =  4i(r)exp  {i(H6  +  kz  -  wt)  }  ,  (A. 3) 

where  l  is  the  azimuthal  harmonic  number,  k  is  the  axial  wavenumber, 
ui  is  the  complex  eigenf requency,  and  ip(r)  is  the  amplitude.  The 
components  of  the  magnetic  and  electric  fields  can  then  be  expressed  as 

B  (r)  =  i(H/r)A  (r)  -  ikA  (r)  , 
r  z  v 

B  (r)  =  ikA  (r)  -  (3/3r)A  (r)  ,  (A. 4) 

r  z 

VO  =  (l/r)O/3r)(rA0)  -  (i/r)*Ar(r)  , 
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and 


Er(r)  =  -(3 /3r)$(r)  +  (iu>/c)Ar(r)  , 


EQ(r)  =  -(ii/r)4>(r)  +  (iu»/c)Ae(r)  , 


(A. 5) 


E  (r)  =  -ik4>(r)  +  (Wc)A  (r)  . 
z  z 


Making  use  of  Eq.  (A. 3),  the  Maxwell  equations  (A. 2)  can  be  expressed  as 

.2 


13  r  It  ~  \  +  p“)<Kr)  =  -4itp(r)  , 


2\- 


r  3  r  3  r  2 
r 


(A. 6) 


13  3  l  +  1,2 


+  i  ,  i  \  ;  .  ,  2ii.  2  /  \  4n  ;  ,  . 
2~  +  P  )  Ar(r)  -  -Y  A0(r)  =  -  —  Jr(r)  , 


4tt 


3  r  3  r 


(A. 7) 


/  1  3  3  l  +  1  ,  2  \  ;  ,  ,  .  2il  ;  ,  »  4ti  ;  ,  , 

(  7  I?  r  37 - 2  +  P  )  Vr)  +  ~  Ar(r)  =  -  ~  Je(r)  ’ 


(A. 8) 


and 


/  1  3  3  fT  ^  2%  :  .  .  4*  -  ,  . 

I  V  r  »  -  T  +  P  )  Vr)  '  '  —  Vr>  ' 


(A. 9) 


where  p  is  defined  by 


2  2.2.2 

P  =  u  /c  -  k 


(A. 10) 


Introducing  the  new  potential  variables 


A+(r)  -  Ar(r)  1  iA0(r) 


(A. 11) 


in  Eqs.  (A. 7)  and  (A. 8)  gives 


13.  3_ 

r  3r  r  3 r 


U  *  1)  2) ;  _  4n  ; 


+  P 


A+(r)  =  -  ~  J+0r)  , 


(A. 12) 


where  the  current  densities  J  (r)  are  defined  by  3.(r)  -  J  (r)  +  iJ  (r) 

i  z  T  v 

In  many  charged  particle  beam  applications,  the  perturbed 
charge  density  and  axial  current  are  related  by 


\  ~ 
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Ja(r)  -  ncp(r)  ,  (A. 13) 

where  n  is  a  constant.  We  therefore  restrict  the  subsequent  analysis 
to  the  case  where  Eq.  (A. 13)  is  satisfied.  After  a  careful  examination 
of  Eq.  (A. 6),  (A. 9),  and  (A. 13) ,  we  find  that  the  axial  component  of 
vector  potential  is  linearly  proportional  to  $(r),  i.e., 

Az(r)  =  n$(r)  .  (A. 14) 

Making  use  of  Eq s.  (A.l),  (A. 5),  and  (A. 14),  the  appropriate 
boundary  conditions  at  r  =  R^  are  given  by 

i(R  )  -  A  <R  )  =  A  (R  )  =  [(3/3r)(rA,)l  _  _  =  0  ,  (A. 15) 

c  z  c  d  c  r  r  -  n 

c 

provided 

kc  -  nco  ^  0  .  (A.  16) 

Here,  Rc  is  the  radius  of  the  conducting  wall. 

In  a  vacuum  waveguide,  where  p  =  0  =  the  solutions  to 

Eqs.  (A. 6),  (A. 9),  and  (A. 12)  are  given  by 

4>(r)  =  (l/n)Az(r)  =  SJ^Cpr)  ,  (A. 17) 

and 

A±(r)  =  a+JJl±1(pr)  ,  (A. 18) 

where  ♦  and  a+  are  constants,  and  J^(x)  is  the  Bessel  function  of  the 

first  kind  of  order  l.  After  some  straightforward  algebra  that  makes 

use  of  the  definitions  A  «  (A  +  A  )/2  and  A.  =  i(A  -  A,)/2, 

r  -  +  d  -*  + 

and  the  boundary  condition  A  (R  )  =  [(3/9r)(rA  )]  ■  0  in  Eq.  (A. 15), 

we  r  k 

c 

we  obtain  a+  =  a  in  Eq.  (A. 18).  Therefore,  the  radial  and 
azimuthal  components  of  the  vector  potential  can  be  expressed  as 


V* 
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Ar(r) 


Vr> 


-  a+iJll(pr)/pT  , 


ia+J j (pr)  , 


(A. 19) 


for  kc  -  nu>  4  0.  In  Eq.  (A.18),  the  prime  (')  denotes  (d/dx)J^(x). 

For  an  electromagnetic  wave  described  by  Eq .  (A. 19)  with  a+  j*  0, 
the  boundary  condition  A^CR^)  =  0  implies 


uj2  ,2  “in 

~i~ k  -  ~r 


(A. 20) 


where  a  is  the  nth  root  of  =  0*  Equation  (A. 20)  obviously 

satisfies  Eq.  (A. 16).  Moreover,  from  the  boundary  condition 

( R  )  =  0,  we  also  note  that  the  amplitude  $  satisfies  $  =  0, 
c 

which  is  consistent  with  the  definition  of  the  transverse  electric  (TE) 
mode  polarization  characterized  by  Ez(r)  *  0. 

Similarly,  for  4>  4  0,  we  find  the  transverse  magnetic  (TM)  mode 
dispersion  relation 


^  -  k2  =  -^  ,  (A. 21) 

c  R 

c 

where  B„  is  the  nth  root  of  the  Bessel  function  J„(B„  )  =  0. 

Jin  i  £n 

Since  t  0,  the  constant  a+  =  0  is  required  to  satisfy  the 

boundary  condition  Ag(Rc)  =  0,  which  is  consistent  with  Bz(r)  =  0. 

From  the  Lorentz  gauge  condition  in  Eq.  (A.l),  we  also  find  n  =  w/kc, 

2  2  2 

and  Eq.  (A. 16)  can  be  expressed  as  k  c  -  w  j  0.  To  summarize, 
for  a  vacuum  waveguide,  the  TM  or  TE  modes  are  exclusively  described 
by  the  potentials  in  Eq.  (A. 17)  or  Eq.  (A. 19),  respectively. 

As  a  second  example,  we  calculate  the  electromagnetic  dielectric 
function  for  the  case  where  the  radial  and  azimuthal  current  densities 


are  given  by 
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Jr(r)  =  -iJQ (r) 


c  6(r  -  V 

Alt  R- 


K  , 


(A. 22) 


where  Rq  and  K  are  constants.  In  this  case,  Eqs.  (A. 7)  and  (A. 8) 
can  be  expressed  as 
,2, 


and 


/  1  8  9  i‘+l  2  \  *  ,  s  211  :  ,  %  6(r~V  „ 

(  7  37  r  97  '  “1"  +  p  )  Ar  r  ‘  ~  Vr)  =  —n K  * 

r  r  0 

•><„  .  «(r-0 


/ 1 

9 _  a _ 

l+l  ,  2 \ 

l  7 

9 r  r  8r 

2  +  "  ) 

(A. 23) 


(A. 24) 


Substituting  Eq.  (A. 11)  into  Eqs.  (A. 23)  and  (A. 24),  it  is  found 
that  the  solutions  to  the  coupled  differential  equations  (A. 23)  and 
(A. 24)  are  given  by 


A  (x)  =  2 


/a±J£±l(pr)  +  b±Nui(pr)  *  R0  *  r  -  Rc  * 

lC±Jn±l(pr)  »  0  i  r  <  R0  • 


(A. 25) 


where  the  constants  a+,  b+  and  are  determined  by  the  boundary 

conditions  that  the  potentials  A  (r)  and  A  (r)  are  continuous  at  r  =  R„, 

t)  r  u 

and  that  A  (r)  and  (9/9r)(rA  )  vanish  at  r  =  R  .  According  to  Eq.  (A. 11), 
ore 

Eq.  (A. 25)  can  also  be  expressed  as 

a-Vi(pr)  +  b-Vi(pr)  +  vWpr)  +  b+Nm(pr)  * 


Ar(r)  = 


R0  "  r  -  Rc  * 

C-J£-l(pr)  +  C+J£+l(pr)  ’  °  -  r  <  R0  * 


(A. 26) 


and 


a_J£_l^Pr)  +  b_N£_i(Pr>  "  a+Jfc+l(pr)  "  b+N£+l(Pr)  » 


AQ(r)  =  i 


R0  <  r  -  Rc  ’ 

c— J£_l (pr)  -  c+J4+i<pr)  »  0  *  r  <  Rq  • 


(A. 27) 


Multiplying  Eqs.  (A. 23)  and  (A. 24)  by  r  and  integrating  from  Rq(I-e) 
to  Rq(1+e),  with  e  -»  0+,  we  obtain  after  some  straightforward  algebra, 

D*(w,k)C+J£+1(0  =  \  K  ,  (A. 28) 


\  - 


\ 


and 
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D’Oa.iOC^J^U)  -  2  K  * 


(A. 29) 


where  the  dielectric  functions  D£(u>,k)  and  D£(w,k)  are  defined  by 


D£(w,k) 


cjt<oj[(o 
Via)*W°  ’ 


(A. 30) 


and 


D~(oJ,k) 


2J  U)J’(ON  ,(£)  -  [J£U)N^(0  +  Jj(C)N£(C)]J£_1(^) 

(A. 31) 


Here,  ^  =  pRQ  and  c,  =  pRc> 

In  a  similar  manner,  we  also  obtain  Eqs.  (A. 28)  and  (A. 29) 
in  circumstances  where  the  radial  and  azimuthal  current  densities 
are  related  by 

c  6(r-Rn> 

J(r)  =  U0(r)  =  -  r“ 


K  . 


(A. 32) 


In  this  case,  however,  the  dielectric  functions  are  defined  by 


2j£(OJ;(c)/*J£+1(0  _ ____ 

D£(m,k)  =  -JjTOjyajH, +1(0  -  7>MON£U)  +  ^(ONf(0)J£+I(0 

1  1  (A. 33) 


and 


u£(oj;(o 


D£(m,k)  =  -  j  (OJ^jU) 


(A. 34) 
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Fig.  1 


Fig.  2 


Fig.  3 


FIGURE  CAPTIONS 

Plots  of  (a)  the  ratio  Rq/Rc  satisfying  Eq.  (86),  and  (b)  the 

E  M 

corresponding  coupling  coefficients  Q  and  Q  for  several 

JC  5  k  S 

values  of  azimuthal  and  radial  mode  numbers,  l  and  s. 

Plots  of  normalized  growth  rate  0^  versus  (k  +  nk^J/kg 
[Eqs.  (82)  and  (83)]  for  U,s)  =  (1,3),  vb  =  10.  v/rb  =  0.02, 

A  ^  =  0.01,  A/y^mc  =  0.002,  and  Rq/rc  satisfying  Eq.  (86). 

Plots  of  normalized  TE  mode  growth  rate  0^  versus  (k  +  nkQ)/kg 
[Eq.  (82)]  for  rq/rc  =  i/a£+i  s>  ^  =  1  and  several 

values  of  s,  (b)  s  =  3  and  several  values  of  l,  and  parameters 
otherwise  identical  to  Fig.  2. 
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Fig.  3(b) 


